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Abstract 

We study the large- iV volume reduction of QCD with adjoint quarks regularized on the lattice. 
Specifically, we use Wilson fermions, and while our d-dimensional lattice has (d — 1) infinite di- 
mensions, the remaining direction is reduced to a point. We perform a weak-coupling one-loop 
calculation of the free energy as a function of the holonomy in the reduced direction, and map the 
regimes in the bare lattice parameter space where the holonomy averages to zero and a Z^y-center 
symmetric configuration is the ground state. For d = 4 and Nf = 1/2, 1 and 2 Dirac flavors we 
see that the center symmetry is intact in a generous regime of the phase diagram that includes the 
chiral point. Thus we see that large- N volume independence of lattice QCD with adjoint Wilson 
quarks works at weak coupling. Interestingly, we find that this is true even if the quark mass is 
quite large, and this opens a path to study the volume reduced large- iV pure gauge theory. Finally, 
we analyze in detail the UV sensitivity of the one-loop potential and show that treating the reduced 
theory as a (d — l)-dimensional effective field theory requires the introduction of certain relevant 
operators that are a subset of those suggested by Unsal and Yaffe to stabilize the center symmetry. 
This means that different regularizations of the volume-reduced theory can be compared only if 
one includes these terms in the action. 



I. INTRODUCTION 



The phenomenon of volume independence of SU(N) gauge theories at large- N was first 
suggested in the influential paper of Eguchi and Kawai lj. It is an example of an orbifold 
equivalence between two gauge theories that are related to each other by projecting out a 
sub-set of fields. In the Eguchi-Kawai (EK) case, one of the theories (let us denote it by 'A') 
is the gauge theory defined at infinite volume, and the other (theory 'B') is the gauge theory 
with one dimension reduced to a point. Clearly, one can obtain B from A by projecting out 
from A all the fields that are not invariant to translations in the reduced direction. Theory 
A can also be obtained from B with a less transparent projection that involves a Zm global 
symmetry whose order parameter is the holonomy in the reduced direction (the Polyakov 
loop). There is a necessary general condition that is required for orbifold projections to 
become large-A equivalences: the symmetries defining the projection must be symmetries 
of the ground states in the corresponding theories Thus, for volume independence to 
hold, translation symmetry and Zn symmetry need to be intact. This was indeed shown 
already in [l| (technically this was done using Dyson Schwinger equations of Wilson loops). 
This can also be seen in the approach of Ref. 3] and, to our knowledge, that reference is 
also the first to emphasize the role of translation symmetry. 

From the lattice perspective volume independence can have an obvious importance; it 
provides a potential way to significantly save computational effort in numerically solving the 
arge-iV limit of certain gauge theories. Unfortunately, however, and as was discovered in 
4], l5|, in the physical case of QCD (or for that matter of pure Yang-Mills) the Zn symmetric 
vacuum is unstable, and the center symmetry is spontaneously broken in the continuum 
limit of the volume reduced theory. In retrospect this is not surprising since QCD has a 
deconfining transitions when its temperature is large enough, or equivalently when its four 
volume is small enough. Therefore shrinking the volume below some critical value leads to a 
Zn deconfining- like transition jfj. Several ways to circumvent this problem were suggested 
over the years and for a review on these we refer to 7| and to 8| . For an updated state 
of affairs we refer to the review section of the first reference in 9fl, and for an intuitive 

n 

explanation we refer to the second reference of |9l. 

n 

From here on we focus on the way Ref. [2J suggested to overcome the EK instability: 
adding to QCD light adjoint fermions that have periodic boundary conditions in the com- 



1 



pactified direction. We shall not discuss at length the relation of this theory to physical 
QCD, but do wish to mention that it is interesting phenomenologically. Its large- N limit is 
expected to be equivalent to the large- N limit of QCD with fermions in the antisymmetric 
representation. The latter theory, in turn, is a generalization of physical, 3-color, QCD to 
a large- N limit in which fermions are truly dynamical and back-react on the gauge fields. 
Thus, it is a complementary limit to the standard large- N 't Hooft limit, where the fermions 
are in the fundamental representation and are sub-leading dynamically. For more details we 
refer to Ref. 



lOj and to the short summary in the forthcoming publication 



12 . 



What the authors of Ref. [2] showed is that in the presence of massless adjoint quarks, 
the Ztv symmetry remains intact down to zero volume. This was done by calculating, in 
weak coupling, the one-loop potential V as a function of the eigenvalues a =i,2,...,JV of the 
holonomy in the small direction. Minimizing V({8 a }) they saw that a invariant vacuum, 
with a uniform eigenvalue density, is the ground state. Note that Ref. [2] was formulated 
in the continuum, and due to asymptotic freedom, the calculation is reliable because the 
length of the compactified direction is small compared to the dynamically generated scale 
in the theory. 



A related work to Ref. \M is that of Ref. 



121 



ll| . There, the authors considered a 3D theory 



describing a unitary adjoint scalar field coupled to gauge fields and to adjoint fermions. 
When projected by a certain center symmetry, the action of this theory was shown to give the 
action of four dimensional QCD with adjoint fermions defined with one compact discretized 
direction. All the other spatial directions were kept in the continuum. The one-loop potential 
of this theory was calculated in perturbation theory, and dimensional regularization with 
minimal subtraction was used to define certain linearly diverging integrals. The result was 



surprising and the ground state was found to generically break the Zm symmetry - in contrast 
to what one might naively expect relying on the analysis of Ref. [2j. 

The result obtained in Ref. [U[ prompted us to perform the study we present in this 
paper and to ask three simple and related questions. First we ask how can one understand 
the apparent difference between the results of Ref. U| and those of Ref. [2]. Second, we ask 
whether both these studies easily carry to a lattice regularization of the field theory where 
all directions are discretized. More precisely, we ask if it is actually possible to reduce one of 
the lattice directions to have a single site, while leaving all other to have an infinite number 
of sites. This is different from the construction of Ref. in which the compact direction 
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was a continuum (i.e. had an infinite number of sites and a zero lattice spacing). This is 

n 

also different from Ref. which had a single site in the compact direction, but where the 
spatial, three-dimensional, space was a continuum. Third, we wish to know what are the 
bare lattice parameters for which we can expect such volume reduction to be valid. This 
is a question of practical importance since if the adjoint quarks need to be very light for 
reduction to work, then it may be impossible to study them in Monte-Carlo simulations. On 
the contrary, if the vacuum is symmetric for heavy fermions as well, then the suggestion 
of Ref. 2I can be considered as a way to probe pure Yang-Mills at large- N. 

To answer these questions we generalize the calculation in [2[ to the lattice. This means 
we need to pick a fermion discretization and our choice is to work with Wilson fermions. 
The prime reason is that this is the fermion discretization that we use to study this theory 
with Monte-Carlo simulations in the companion publication fl2 |. Our calculation can also 
be considered as a generalization of the pioneering Refs. % |5[ in two ways. Firstly, we 
generalize them to the case where only one of the euclidean directions is reduced to a point, 
and secondly, we add dynamical adjoint fermions to the theory. (The first time these works 
were generalized to a singly reduced direction was already done in but since this was 
not published, we repeat some of its steps here). 

The following is the outline of the paper. In Section [III we define the action of the theory 
that we study, and prepare for the one-loop perturbative calculation by fixing axial gauge. 
In Section [Till we calculate the contribution of the gauge fields to the one- loop potential. As 
mentioned above, this section is very similar to what already appears in [131 ] . In Section HVl 
we calculate the contribution of the adjoint fermions to the potential. The form of the 
resulting potential is summarized in Section [V] where we also remark on several of its 
properties and on the way its parameters should be chosen in order to get continuum- like 
equations like the ones that appear in Ref. Jjj. Before we minimize the potential we pause 
in Section [VI] to discuss the way the one- loop potential was calculated in Ref. H[ , explain 
the meaning of the results obtained in that work and why they do not signal any problems 
for large- N volume reduction. This is tied with the way Ref. [U] treats the UV divergences 
of the one- loop potential, and we therefore analyze these in detail. In Section [VIII we return 
to our lattice regularization and present a map of its phase diagram along the quark mass 
axis for the case of an isotropic lattice. This is generalized to the case where the lattice 
spacing in the reduced direction a t is different from the lattice spacing in the unreduced 
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dimensions a s . In that case we map the phase diagram in the plane of the quark mass and 
the ratio a s /a t . We summarize our findings in Section HXl 



II. THE ACTION AND AXIAL GAUGE FIXING 

The starting point of this section is the action of ci-dimensional QCD with Nf Dirac 
fermions regularized on a <i-dimensional hyper-cubic lattice. Out of the d space-time dimen- 
sions, d — 1 have an infinite number of sites and one has a single site. The gauge action is 
the standard Wilson unimproved action, and the fermions are chosen to be of the Wilson 
type. Also, we allow for anisotropic lattice spacings and so have two gauge couplings and 
two hopping parameters. Those that correspond to spatial plaquettes and hopping terms 
carry the super /subscript "s" while those corresponding to time- like quantities will have 
the super /subscript "t". In pri ncip le, one can add an anisotropy in the Wilson coefficients 
(the r), but this is not needed [14J. The boundary conditions on the compact direction are 
periodic for all fields. The action is then 

A ^gauge ^gauge ^Fermions (2-1) 

2N 

^gauge = ~T~ R- e £ Ux,iU x+ ijUl +j<i Ulj, (2.2) 
A s x 

i<i£[l,d-l] 

2N 

4auge = ^-Re J2 ^i'r,r,.,J J ..,rL. (2.3) 

t X 

i£[l,d-l] 

^Fermions = V> D W $ , (2.4) 

£(i- 7.) u£i s y , x+ i + (i + 7.) u% ^ 

i 

-7o)^ G ^ + (l + 7oX?^]. (2.5) 

Here the (d — l)-dimensional site index is x, and the indices i, j denote the d — 1 transverse 
coordinates. We denote the compact direction by 0. Also, while U X>IM are the usual link 
matrices, then denote their adjoint representation. Lastly, we denote the 't Hooft 
couplings gl s N by A t>s , and the 'hopping' parameters by « S)t . 

We begin by fixing a gauge in which all the temporal links are diagonal. This is done by 
writing each temporal link as a unitary conjugation of an x-dependent diagonal matrix that 
we denote by e*** = diag (e ilpi , e iip2 , e itpN ) x : 

U xfi = V x e l ^Vl (2.6) 



If we keep V x G U(N) and (f a > ip b for a > b, then this is almost a one to one mapping 
because given a pair of matrices U Xj o and e l(Px , there is more than one matrix V x that satisfies 
Eq. (12. 6ft . Specifically, any two matrices V x ^ and that are related to each other by the 
right multiplication 

V® = V® x A,, (2.7) 

with A x an arbitrary unitary diagonal matrix, will give, for a given set of eigenvalues {(fi a }, 
the same matrix U x $. This JT {7(1)^ symmetry is a residual gauge symmetry which we will 
not fix, but do not find it to be important at large- iV (see below and Pij]). Using Eq. ( 12.61) . 
the measure of the temporal links on each spatial site becomes 



DU X>0 = / DV X 

U{N) Ju(N) 
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, 2tt 



n-w^) (2.8) 

a<b ^ ' I v a> v b 



The next step is to use the U(N) spatial gauge symmetry of the path integral and perform 
the following change of variables: 

U Xti -> VjU x>i V x+h (2.9) 
^ -> (2-10) 
4^^K tG (2-H) 

This makes the path integral and any gauge invariant operators independent of the diago- 
nalizing matrices V x , and we are left with the following gauge- fixed path integral 

x,i ^ U ( N ) a J a<b ^ ' 

(2.12) 

This will be our starting point for the one-loop calculation. 

III. THE GLUONIC CONTRIBUTION TO THE ONE LOOP POTENTIAL 

In this section we calculate the gluonic contribution to the one loop potential. This is, of 
course, not the first time it is done, and some related references are j^],!^, Q . We begin by 
noting that the vacua 

u£ t = (3-1) 

</% = (3-2) 



are classical maxima of the action that are degenerate for all choices of 9. In this section we 
integrate the Gaussian fluctuations around these vacua to break this degeneracy. We do so 
by writing 

jj* = gat + ^ _ 1 + 0(A 3 ); (3 3) 

= e a + r x - (3.4) 

After some familiar algebra one finds that the second order contribution to the spatial 
plaquette terms in the action, A s gange , becomes 

O 7V" 

/is,2nd-order _ \ ^ BpTr \ A A . . i A A . i A A . . i A . A . 

gauge \ / iLcil | n i,i' 1 l+ij ~ - rl -x,i- r '-x,] ~ ■ n -x+],i- rl x+ij ~ ■ rL x+j,%- rL x,% 



Xs 

X,l<J 

_A A A ..A _1A 2 --A 2 ~-A 2 --A 2 



(3.5) 



and that the temporal part of the action, ^4g aug e, takes the form 

Efei-« 2+4 Ei^i 2sin2 ( 



4<,2nd-order _ ^ 
gauge — ^ 



(3.6) 



Let us make the following three remarks: 



We have not taken into account the way the measure of the spatial link matrices 
depends on A and so have assumed it to be flat. The leading correction is in fact 
quadratic in A and will formally contribute to the one loop potential [3| . Nonetheless, 



because we are at weak coupling, the quadratic pieces that we do include in Eqs. (13. 5j) - 
(13.61) are large by factors of l/Xt,s ^> 1 compared to those quadratic pieces that we 
neglect. Thus, truncating the action as we do in Eqs. (l3.5H3.6p is indeed a consistent 
approximation. 

At the quadratic level the <fi fields and the A fields decouple. Because the action of 
the former does not depend on 6, we ignore it from here on. 

Fluctuations along the diagonal components of A are completely flat. This reflects the 
residual gauge symmetry discussed above (see Eq. (12.71) ). The action of these fields 
will become non-flat once we appropriately choose a gauge, but because there are only 
N such fluctuations, integrating over them will make only an O(N) contribution to the 
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effective potential. Since we are mostly focused on the large- N limit of the potential, 
such a contribution will be 0(1/ N) suppressed compared to that of the off-diagonal 
components of the gluons and as a result we can ignore these fields. 

We are now left with integrating only over the off-diagonal components of the spatial 
gauge fields A a jf . To do so we Fourier transform according to 



A% = 1 (f) e^PAt, (3.7) 



d-l 

ikx+iki/2 a 

\ki\<TT V Z7r / 

and find the following form of A^- order = ^nd- order + A *,2nd- order 

n r fdkY' 1 



*S£«~ = / (£) " EE [A&M$(k)A%] , (3.8) 

s ij a^b 

Mij(k) = 5 l3 (y2sm 2 (h/2) - sin 2 (^/2) + ^ sin 2 



-(1-dij) sin(^/2) sin^). (3.9) 

Since A c ^ i = (A b ^ k /)* there are only (d — 1) x iV(iV — l)/2 independent coordinates for each 
value of k and we write 

s j \ / ^ a<b 

(here we used M^(k) = Mf b {— k) = M i 6 j l (/c)). Integrating over A gives the following contri- 
bution to the one-loop potential (the second term comes from the measure of the 9 fields - 
see Eq. (TXED ): 

w*) = E / (tY 1 log K et <H ~ E lo s ( sin2 (^y^)) • ( 3 - n ) 



With simple linear algebra one can show that the (d— l)-dimensional matrix (M ab {k)) 
has (d — 2) identical eigenvalues equal to £V sin 2 (A;j/2) + ^ sin 2 ^^_=^ J an d one eigenvalue 
equal to ^ sin 2 ( ^ a - 6b j sj]. 1 Dropping ^-independent terms, it is easy to see that the 



1 To sec this write (M ab (k)).. = (jT Jl sin 2 (fc l /2) + sin 2 (^^)) % - (A/ ab (fc)) and diagonalize M(k). 

The latter has the structure My = s^s.,- with = sin(fcj/2) and so has one eigenvector proportional to vq — 
(si, S2, ■ ■ ■ ! Sd-i), with a corresponding eigenvalue equal to J2i s f- There are also (d— 2) eigenvectors that 
are orthogonal to vq and that therefore correspond to zero eigenvalues. This gives the set of eigenvalues 
for the matrix M discussed in the text. 



7 



contribution of the latter to Eq. (13. lip cancels the second term in Eq. (13.111) , and as a result 
one obtains the following form 

w«) = tS / (^""'"sfe^^ + r 8 " 1 ^^-)}' (3 - 12) 

By sending the oil loop potential in this fcnn cold have been anticipated in 

advance. Also, the L t > 1 version of Eq. (13.121) already appears in 3]. 2 Eq. (13.121) is indeed 
suggestive that, in the invariant vacuum, 9 a — 9 b is what plays the role of the fourth 
component of the gluon momenta. This is a generic feature of the way single-site reduced 



models embed the first Brillouin Zone of the large volume theory into color space 



16| 



D xy 



IV. THE FERMIONIC CONTRIBUTION TO THE ONE LOOP POTENTIAL 

To add adjoint Wilson fermions we use the form in Eq. ( 12. 5p . At one loop we set U x ,i = 1, 
and U^q = 5 a b e t9a . Then D is diagonalized in color space and we have 3 

5 X y - K s [(1 - Ji)5y )X+ i + (1 + 7i)<W-t]| 

-«*{(! " lo)e^- eb) + (1 + 7o)e- i(ea - eb) } . (4.2) 
In momentum space, D has eigenvalues given by (here we denote 6 a — 6 b by 9 ab ) 

{£[(1-70** + a +70^]} 

-^{(l-7o)e^ + (l + 7 o)e-^}, (4.3) 

which can be written as 
D p = 1 - 2k s ^2 cos 

h + 2K t cos {6 ab ) + 2ik s J2 Mkhi + 2 % «* sin(# afc ) 7o , (4.4) 



Dp — 1 K s 



2 Although with an overall factor of two which seems to be redundant. While unimportant for Ref. [3] 
(that studied only the pure gauge case), this factor is crucial for us, as we are comparing the magnitude 
of Eq. (|3.12[) to the contribution of the fermions to the one-loop potential (see section ITV)) . 

3 Here we work in a basis where the color indices A and B of the adjoint representation link (U G ) AB are 
composite indices that correspond to pairs of fundamental indices: (U G ) AB = (U G ) ab cd — U ac U£ d . This 
means that the classical configuration of the time-like link is given by 

(U G ) x0 \ = S ac S bd e^ a - e ^ ee 6 ac 6 bd e* eab . (4.1) 

' J ab.cd 
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or as 



D p = (l-2« - (d-l)-2« t )+4K a + — sin 2 (9 ab /2)^j +2w s 



V sm(^)7, + ^sin(r b ) 7 o 



(4.5) 

Here we defined S = ^ i sin 2 (/cj/2). In d space-time dimensions, and for Afy Dirac flavors, 
the determinant of D p is given by 

2 



det D r , 



j((l-2«,(d-l)-2Ki)+4«. (^ + ^ sin 2 (0 a6 /2)^ 



+4k 2(5_+(^) sin 2^ 



2 (d/2-i) A r / 



(4.6) 



where we have also defined S2 = J^sm 2 fcj. 

Note that the calculation in this section ignores the fermionic zero modes. As shown 
in Ref. 17J these are important in the supersymmetric case (i.e. when Nf = 1/2 and for 
massless quarks). In that case the one loop potential at the vicinity of the invariant 
phase is zero, and the subleading, O(N), contribution of the fermionic zero modes is what 
is left. Since we are not particularly interested in the SUSY theory, and because away from 
the Zn invariant state, SUSY is explicitly broken, then we can always neglect these zero 
modes. 



V. RE-CAP AND CONNECTING THE LATTICE PARAMETERS TO THE PA- 
RAMETERS OF THE CONTINUUM CALCULATION OF REF. |lj 



To conclude, the following is the 1-loop potential of the full theory: 

V g (k, X s ,X t ) 



a^b 



V f (k,N f ,K s ,K t ) 



V g (k, \ s ,X t ) 

V f (k,K s ,K t ) 



' M Xs . 2 f9 ab 

s+ \~ t sm It 



(d-2)/2 



„ / gab 

(1 - 2K s (d - 1) - 2«t) + Ak s ( S + — sin 2 — 

K s V 2 



2 d/2-i A r 



: 4k 2 |S 2 4 (^) sin 2 ^ 



(5.1) 
(5.2) 



(5.3) 
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i=l ^ ' 

d-1 

5 2 = sin2 (**)■ ( 5 - 5 ) 

i=l 

One way to see that Eq. (15.11) makes sense is to set d = 4 and Nf = 1/2. In that case the 
power of the outermost brackets in both Eq. (15.21) and Eq. (I5.3P is equal to 1 and so if we 
take the large-iV limit in the invariant phase, and appropriately tune the bare lattice 
parameters to have massless quarks in the continuum limit, one can see that V(9) —>■ 0, as 
it should, due to supersymmetry. 

Let us now show how to tune the parameters in Eq. ( 15. ip to formally get continuum- 



like equations like those obtained in Ref. III]. 4 There, the authors work with a model 



which naively looks like a continuum version of the reduced model that we study here. To 



see the connection between our Eq. (15.11) and Eq. (3.14) of Ref. llj we re- introduce the 
spatial and temporal lattice spacings a s and a t . This is done by writing ki = a s pi and 
dividing T> g and T>f by appropriate powers of a s such that S in the numerator and S2 in 
the denominator will both get divided by a 2 s . Then, if we take a s — > but keep a t fixed, we 
have (4S/a%, S^/af) — > J2iPl = P 2 - This procedure turns T> g /Vf into the following form 5 

(d-2)/2 

LJg a t =fixcd Y At<1 « ' 

7 / (l-2/t s (d-l)-2/Cj) 

(5.6) 

[ — 1 

Comparing Eq. (15.61) with Eq. (3.14) of [11], we see that the following relations connect 
the bare lattice quantities of our calculation (i.e. K t , k s , A t , A s ) with the bare quantities of 
continuum-like calculations of the sort of [ll] (i.e. the quark mass in units of the temporal 



+ ^ + Si sin 2 (¥)) +U+ (^) sin 2 6* 



lattice spacing, a t m): 



.6 



at = a s \ T~> (5.7) 
V A. s 



A related calculation for the pure gauge case was performed in Ref. 

5 Here we also divide the content of the square brackets in the denominator by 4k^. This gives a potential 
that differs from the original one by a 9 independent term, which we ignore. 

6 There may be typos in some of the numerical coefficients of the relevant equations of Ref. [Ill ] , and for 
example it seems that the argument of the (sin) 2 term in the equation at the bottom of page 10 should 
be divided by 2, while its coefficient should be multiplied by a factor of 2. We think this leads to an 
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K t 



a t 



2(d - 1) + 2(a t m + Jj*) 



(5.8) 
(5.9) 



Thus, naively, one might expect that as we take A s /A 4 — > we should recover the results 
reported in Ref. There, the vacuum was seen to be a state that breaks the center 

symmetry, except for a small, physically uninteresting, window of very heavy quarks with 
a t m e [0.6, 0.8]. 7 



In the next section we discuss the results of Ref. 
is wrong. 



Ill ] and show that such an expectation 



VI. THE VOLUME-REDUCED THEORY AS AN EFFECTIVE FIELD THEORY 



In this section we explain why the results of Ref. do not signal a problem for large- 
N reduction of the theory. To show this it is sufficient to focus on the gauge field sector 
- later, in Section IVI El we will show how fermions modify the discussion. We therefore 
begin by writing the gluonic contribution of the action studied in Ref. It is the naive 
continuum limit of Eqs. (I2.2l) - (l2.3p taken only in the spatial dimensions. Thus, it describes a 
D-dimensional gauge theory coupled to a non-linear unitary, adjoint sigma field Q G SU (N) 
and its action is 

-, D , 2 D 



S' 



one— site 



dPx tr 



2g 2 



(6.1) 



DpQ. = d^l + i [An , fi] 



(6.2) 



The case studied in Ref. 



ll| had D = d - 1 = 3, g 2 = 2X s /(Na s ), and f 2 = 2N/(X t a s ), but 



for our discussion it is convenient to leave g and / as free parameters. 

Had we allowed for L t sites in the compact discretized direction, then the continuum gauge 
fields A^ix) and the unitary field Q would acquire an additional index t — 1, 2, . . . , L t . In 



extra pre- factor of v8 to the first term in the square brackets in the equation at the bottom of page 11 
and in Eq .(3.14). This should not change the conclusions of 11] since this factor strengthens the Zm 
destabilizing gluon contribution to the one-loop potential. 
7 We are working with symmetric Wilson parameters that are equal to one. 
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that case the field flt(x) transforms under a bi-fundamental representation of the SU(N) 
groups generated by A^(t) and A^{t + 1). The action is then 



S 



it— sites 



E [ dDxtT it F l (Mt)) + y E 

t=l ^ L fi,v=l M=l 



(6.3) 



D M a = <9,A + iA„{t)n t -in t A lx (t + l). 



(6.4) 



An action of the form of Eq. ( 16. 3j) has been studied in the context of a low energy ef- 
fective action for strongly-coupled Higgs models [18[, deconfinement of Polyakov loops in 
four dimensions 13 , and of deconstruction of five- dimensional theories as candidates of the 
standard model 20j. Both former studies worked with D = 3 and L t — 1, while the latter 
studied D = 4 and L t > 1. 

Below we first discuss several features of 5 ,it_sltes that are important in our context 
( Sect ions I VI AITVTCl) . then discuss the results of Ref. [jjj (Section IVIDI) . then analyze the 
way fermions affect our arguments (Section IVI El) , and finally summarize in Section IYl Fl 



A. Non-renormalizability 

As noted by Refs. 2l3], the theory defined by Eq. (\G.2\\ with D = 4 is non-renormalizable. 
For D = 3 this nonlinear sigma model is also not expected to be a renormalizable field theory, 
at least not in perturbation theory (for example see the remarks in the second reference of 
H or in Refs. JjJ)- 

In fact, dropping the contribution of the gauge fields to Eq. (16.21) and setting L t — 1, we 
obtain the two- derivative term in the chiral Lagrangian for an SU(N) x SU(N) — ► SU(N) 



breakdown scheme 



2l| . For L t > 1 the model is a more complicated sigma model describing 



the spontaneous breaking of SU(N) 2Lt to SU(N) Lt 20[]. While the sigma model for L t = 1 
and D = 4 is certainly non-renormalizable (again see 21|), we are not aware of detailed 
studies of this issue for the other cases. A simple analysis, however, shows that for the case 
of interest in this paper (D = 3 and L t = 1) the action in Eq. (16. 2p . when expanded in 
terms of the 'pion' field 7r defined by fl = e ln , will radiatively generate four derivative terms, 
which will then generate eight derivative terms, etc., and that all these terms will be linearly 
UV- divergent. While a similar effect happens at the one- loop level for D = 4, then the lower 
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UV divergences in D = 3 generate the higher derivative terms only at the three-loop level. 
The one-loop potential in our case, however, is UV sensitive already at one-loop (see below). 

In any event, since Eq. ( 16.31) is non-renormalizable, we can understand it as an effective 
field theory (EFT). 

B. UV sensitivity of the one-loop potential for low values of L t 

The one-loop potential of our theory contains UV divergences. Denoting a momentum 
cutoff by A, we should expect a leading A D divergence which can be removed by subtracting 
from the one-loop potential its value at 9 = 0. Importantly, however, there are sub-leading 
divergences that depend on 9, and that can (and do) affect the breaking scheme of the 
symmetry. To show this we first identify the leading divergences in our lattice calculation, 
where A ~ l/a s . This is done for L t = 1 and D = 3 in Section IVl B II We then discuss in 
Section fVI B 21 the general case of Lt > 1 and of D spatial dimensions. 

1. The case of Lt = 1 and D = 3 

The starting point is to note that the dimensionful one-loop potential is given by multi- 
plying V{9) from Eq. (15. lj) by 1/ (a 3 a t ). Next, we subtract from V{9) its value at 9 = and 
obtain 

' W) _ no)r 'rf (*y iog ( 1 + ±!!!lMY (6 . 5) 

To identify the UV divergence we expand the log in 1/p 2 and obtain 

Dropping a ^-independent constant this can be written as 

-i- [V{9) - U(0)]™ divergent ~ 1 |tr fi classical | 2 x 1 (6.7) 
where we re-identify the classical values of the unitary diagonal holonomy field Q as 

( ^classical ) a fc ^ $ab- 

(6.8) 
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Thus we see that the (mass) 2 term of the classical holonomy is linearly UV- divergent at 
one-loop. If we add fermions this will not change, except for one case: for Nf = 1/2, m — 0, 
and in the vicinity of the Zm invariant vacuum, the one-loop potential is identically zero 
due to supersymmetry. But away from this vacuum or for other values of m and Nf, the 
^-dependent linear divergences remain after the introduction of fermions. 

2. The case of general L t and D 

In general, the UV sensitivities of the one-loop potential can be seen from its expansion in 
terms of bubble diagrams. 8 For that we focus on the quartic interaction between the gauge 
fields and the Q fields - see the second term in Eq. (16.21) . For general L t these interaction 
terms are given by 

g 2 ftr (n t A,(t + 1) njiyt)) = g 2 f J>f Aftt + i^A**®, (6.9) 

abed 

which we depict pictorially in Fig. [Tj 




=(gty 



FIG. 1: The vertex of a term of the form of Eq. (TTJ). 

Gauge invariance tells us that the only terms that can be generated radiatively are made 
out of the Polyakov loop P(x) defined by 

P{x)= (jl^t(x)\, (6.10) 
and the center symmetry allows only operators of the form 

|tr P(x)\ 2 , |tr P 2 (x)\ 2 , |tr P 3 (a;)| 2 , . . . (6.11) 
8 We thank C. Hoyos for suggesting to us this diagrammatic way of thinking about the one-loop potential. 
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(we ignore the fact that operators like (tr P) N are also allowed by the center when the gauge 
group is SU (N), because at large- N there should be no difference between SU(N) and U (N) 
and the center of the latter gauge group does not allow such operators). 

These symmetry restrictions can be seen by drawing all the possible disconnected bubble 
diagrams that contribute to V(9). For example, the diagram that gives rise to the first 
operator in Eq. ( 16. 1 If) with L t = 8 is given in Fig. [2j A simple counting show that this 




FIG. 2: The bubble diagram that gives rise to the term |trP| 2 with Lt = 8. The numbers next to 
the vertices denote the t index of the Q and fit fields on the external legs, while the numbers on 
the gluon lines denote the t indices of the A fields that flow in the loop. 



diagram scales like 



± dv jM^ y r^ (6 , 2) 



which is finite in the a s —>■ limit so long as 

D < 2L t . (6.13) 

Ref. [13J already showed this by explicitly calculating the one-loop potential for the pure 
gauge theory as a function of L t . 

Thus, for D = 3 and L t = 1 we get the linear divergence discussed above, while for 
L t > 1 one expects no divergences at all. For D = 4, one expects a quadratic divergence for 
Lt = 1, logarithmic for L t = 2 and no divergences for L t > 2. Indeed, this power counting 



led Ref. [20|] to consider the L t > 3 system as a model for a four dimensional theory that 
dynamically generates an extra dimension, and whose Higgs mass is not UV sensitive. 

Terms in the one-loop potential that correspond to terms in Eq. (16.111) with k > 2 powers 
of P inside the trace (k windings of Polyakov loops) are generated by bubbles diagrams that 
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are similar to the one in Fig. [21 but with kL t vertices that are ordered cyclically. Thus, the 
UV contribution to these diagrams is of O (f ^77) and is finite as long as 

D < 2kL t . (6.14) 

This means that these contributions are all convergent for L t — 1 and D = 3, while for 
Lt — 1 and D = 4, the k = 2 contribution is logarithmically divergent. 

The fact that the 9 dependence of V(6) is sensitive to the UV regulator when L t is small, 
is, at first glance, quite surprising. Does this contradict the result of Ref. where no 
divergences were observed? The answer is no and the reason is that the calculation there 
was defined in the continuum of the compact direction and effectively had L t = oo. Thus 
in this case the spatial integrals we discuss above are convergent. Technically this is due to 
the infinite sum over the Matzubara frequencies corresponding to the compact dimension. 
Denoting the length of this direction by R, this infinite sum yields exponential terms of the 
form e -R ' p ! that suppress the higher spatial momenta of 0(1/ R). 



C. Relevant operators that are missing from the original action in Eq. (16. 3j) 

As mentioned above, the non-renormalizability of the action in Eq. (16.31) tells us that we 
can see it as an EFT. This means that we can add to the action all possible operators that 
are allowed by symmetries. The coefficients of these operators are a priori arbitrary and can 
be chosen at will. As the discussion above shows, we can consider terms of the form 

5S W = Jd D x^B 1 \tiP(x)\ 2 + B 2 \tiP 2 (x)\ 2 + ...Y (6.15) 

with coefficient B 12i ... of mass dimension D. Other operators we can add have coefficients 
with negative mass dimensions, for example 

SS {2) = J2 [d D x {d (tr |ZW)| 2 ) 2 + C 2 (tr \D^(t)\ 2 f + ...}. (6.16) 

While we can choose some of the coefficients of these operators to zero, it is important to 
understand that at sufficiently high loop order, they may be generated with UV diverging 
coefficients. In that case, from the EFT point of view, we will need to add them to the 
action as counter terms whose bare couplings contain UV divergent pieces that cancel the 
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divergences, but also finite pieces that will become new low energy coefficients of the EFT. 
This is standard in EFT: going to higher loop order means one is sensitive to more low 
energy constants. 

In our case we are working in one- loop, and so need to add the terms to the action that 
will be generated with UV-divergences at that order. For D = 3 and L t = 1 this means we 
need to add the following term 



The coefficient c\ in Eq. (16.171) will cancel the linear divergence of Eq. (16. 7p . but the finite 
low energy coefficient b\ can be chosen at will. Different choices of b\ correspond to different 
parts of the EFT's space of parameters. Therefore, we see that at one-loop the EFT in 
Eq. (16.21) is defined not only by g and /, but also by b\. Clearly, for the purpose of large- 
N reduction, a preferred choice of b\ will be one that leads to a one-loop potential whose 
ground state is i^Ar-invariant. 

For D dimensions and L t time slices we will need to add terms of the form 



where |_^J is the integer part of x, and when D = 2L t k there is a logarithmic divergence (for 
brevity of notations, we assume that each term in Eq. (16. 18ft is radiatively generated only 
with a leading divergence. This need not be the case in general). 

Before we proceed let us emphasize the following two issues: 

1. Terms of the form of Eq. C I6. 181) need to be added to the action only if one wishes to 
treat the reduced model .D-dimensional EFT. 

But this is not our purpose. What we wish to do in this paper is to check whether the 
theory, as defined in Section [Til an d for a given set of bare lattice parameters X Sjt and 
K Stt , has a vacuum that is symmetric in weak-coupling. We do not need to think 
about our theory as an EFT, nor do we need to take its a s — > limit and worry about 
canceling divergences. If we find that the ground state of our L t = 1 reduced model, 
defined with fixed cutoffs a s and a t , is Zn symmetric, then large- N reduction tells 
us that it is large- N equivalent to the L t = oo four dimensional theory defined with 




(6.17) 




k=l 



(6.18) 
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the same field content, the same regularization, and the same cutoffs. To remove the 
four-dimensional cutoffs we will then tune the lattice parameters according to their 
four dimensional RG flow. According to the large- N equivalence paradigm, this should 
be done only after taking the large- N limit. Thus in this approach (which, in fact, is 
the standard one used to justify large- N reduction non-perturbatively 0, [26]) we 
do not add any counter terms to the action. 9 



In contrast, in Ref. lfj, the authors do treat the reduced theory as an EFT. Specifi- 
cally, minimally subtracted dimensional regularization was chosen to get rid of the UV 
divergences. For the motivation of treating the theory in this way see the introduction 
of that paper. 

2. It is interesting that the terms in Eq. f)6.18p . that we need to add to the action from 
the point of view of EFT, are a subset of the terms that Unsal and Yaffe suggested to 



,8 



As shown there, 



add to the EK model in order to stabilize its invariant vacuum 
when projected back to 4D, these terms change the gauge theory in a way which is 
only sub-leading at large-A. 



D. Connecting the results obtained in Ref. 11] to other regularization schemes 



The lesson of the previous sections is that if we view the reduced theory as a D- 
dimensional EFT, then the action describing it is not given by Eq. f)6.2p ; the latter has 
certain terms missing. In particular, for D = 3 and L t = 1, the most general EFT is given 
by 

d e2 D 



+ (c 1 A + 6 1 )|trP(a;)| 2 !.. (6.19) 



Here we need to tune c\ in a regulator-dependent way, and to choose &i as we wish - different 
choices correspond to different low energy constants of the EFT. 



9 There is of course another reason to add terms of the form |trP fe | 2 with k > 1; such terms can get rid of 
unwanted center symmetry breakdown in regimes of the lattice parameter space where it surely happens 
(see details in Ref. Q and below). In that case, however, they are not considered as counter terms: 
their coefficients need not be tuned like 0(l/a s ) and the choice of the power k is not dictated by the 
considerations in Section IVl Bl 
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Since different regularization schemes subtract the UV divergences in different ways, 
then one can compare regularization schemes only if one adds Eq. (16.171) to the action. For 
example, the regulator used in Ref. JUj was minimally subtracted dimensional regularization 
(MSDR). This regulator sets power law divergences, like the linear divergences in the one- 
loop potential, to zero, and essentially replaces the UV divergence of V(6) by finite 9- 
dependent functions. This fact makes MSDR a 'dangerous' regularization scheme in our 
context - it automatically subtracts the infinity, and in contrast to other regularization 
schemes, it does not make the necessity of the counter terms in Eq. (16.191) manifest. In fact, 
this was already noted in Ref. [23]. 10 

Indeed, because MSDR sets the linear infinity to zero, then there is nothing to subtract 
and the infinite piece of the counter term (ci in Eq. ( 16.191) ) is fixed to zero in this regular- 
ization. Since, however, the last term of Eq. (16.191) was not considered in Ref. fll|. then b\ 
was implicitly set to be zero as well. Consequently, the resulting V{9) was minimized in a 
subspace of the EFT full parameter space, and this subspace is not special in any sense (the 
point b\ = is not protected by any symmetry). Therefore the fact that V(8) was found to 
break the center symmetry in a physically relevant region of the EFT parameter space is a 
result that is particular for the choice b\ = and may certainly change once one explores 
other choices for b\ which correspond to other choices of regulators. 11 

For example, we can choose any 3D lattice regulator for Eq. (16.191) (such as the standard 
Wilson action, or an 'improved' one), or variants of dimensional regularization like power 
divergence scheme j^. If one does not add the counter terms to the actions of these 
regulators, then each regulator makes its own implicit choice for b\. As an example, let us 
choose the regulator to be a lattice and use the standard Wilson action. This is in fact the 
same action that we used in Section [III (While in Section IVI CI we emphasized that our 
lattice reduced model is not taken to be a regulated EFT, we can momentarily depart from 
this point of view, and use it as one). What are the values of c\ and b\ that our lattice 
calculation chooses when viewed as an EFT? In general, if we denote by c[ the coefficient 
of the l/a s term that multiplies |trP| 2 in the lattice result for V(8), then we need to set 



10 



See Appendix C_Lhere. We thank M. Unsal for bringing this reference to our attention. 



11 The action Ref. ll| studied also contained fermions and, compared to Eq. (|6.19p . it has more terms in its 
action (which require an one additional counter term). In the next subsection we discuss this issue, but 
for the arguments in the current subsection it is not essential. 
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Ci from Eq. (16.1 7p to be —c[. Next, since the action that we use in our lattice calculation 
is only Eqs. fl2.2p - fl2.3l) . then had we treated it as an EFT means that we actually chose 
bi = -ci/a s = +c' 1 /a s . 

Our message in this section is that the absence of counter terms leads the regulator to 
implicitly choose different values of the low energy constant b\. This choice is determined by 
the details of the regulator and so even if we fix the physical parameters / and g, then there 
is no reason why the two regularizations will yield the same physical result. In particular, 
it is quite possible that while one regularization sees a invariant vacuum, then the other 
concludes that the symmetry is intact, even if / and g are the same in both. To get 
identical physical results in any two regularizations, one will need to explcitly add the counter 
term to their action, and tune the values of b\ in both regulators in an appropriate manner. 



E. Effect of adjoint fermions on the UV sensitivity of the one-lop potential 

In this section we show how the presence of the adjoint fermions modifies the discussion 
above. The modifications are two- fold. First, because of the fermion propagators are of 
0(1 /p), the one- loop divergences can be of a higher degree. Second, because the Dirac 
operator of the fermions is a first derivative in the compact direction, which becomes ~ 
sin(6 ab ) in the reduced model, then the theory will radiatively generate operators of the 
form |trP 2 | 2 as well as |trP| 2 . 

These facts can be seen in two ways. From arguments of the sort of Section I VI B 1 1 we 
see that the UV sensitive ^-dependent pieces of the fermionic contribution to the one-loop 
potentials are given by 

r b \ . 2 (e ab \ . 2 fQab s] f d 3 P 

— j , or sin 2 i^—j , or sin 2 (8 ab )^ x J — , 

(6.20) 

and a simple rearrangement of the color indices tells us that both |trP 2 | 2 and |trP| 2 are 
generated with a linearly diverging coefficient. This means that the EFT needs to also 
contain the |trP 2 | 2 operators, and is thus defined not only by the quark mass and by /, g 
and bi, but also by the values of a new low energy constant b 2 that comes from a counter 
term of the form 

fermions = f fa+A | tr P 2 (;r)| 2 . (6.21) 



^|v(e)-v(o)|g^z«~E 



3'„ 

06 



a^at 
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Arguments relying on the structure of bubble diagrams, like those of Section IVI B 21 can 
also be used. The bubble diagrams are now generated by replacing the vertex in Fig. [T] with 
a corresponding vertex that connect two quarks and two Qt matrices, and by replacing the 
gluonic loop of Fig. [2] by a quark loop. Since the trace over the Dirac gamma matrices will 
null all diagrams with an odd number of vertices we find that for general D and L t one 
needs to generalize Eq. (16.181) to a sum over k = 2,4,6, ... , fc max where k max = if |_^J 
is even and fc max = ([^J - 1) if is odd. 



F. Summary 

A simple and important result of the previous subsections is the following. If one wishes to 
treat the reduced model as a three-dimensional EFT, then V(9) calculated for the action in 
Eq. ( 12. ip (if we regularize the EFT with a 3D lattice) or for the action Eq. (\G.2\\ generalized 
to include fermions (if we choose MSDR), is missing the following terms 

hissing ~ 6i|tr P\ 2 + b 2 \tiP 2 \ 2 . (6.22) 

Different values of b± and 62 correspond to different points in the parameter space of the 
EFT. Thus, to show that for given values of physical parameters like a t m and Nf, the EFT 
spontaneously breaks the center symmetry, we need to verify that there is no combination 
of bi and 62 that can make the ground state of V(9) center invariant. Since this procedure 



was not included in the analysis of Ref. 11] . and instead, the regulator used there implicitly 
fixed bi = 62 = 0, then it is certainly possible that for the same values of m and Nf there 
is a different point in the plane spanned by b\ and 62 for which V(9) has a Zn invariant 
ground state. 

Let us show that this is very plausible. First, note that because different regulators 
subtract the UV divergences in a way that differs by finite pieces, then two regularization 
schemes will give the same physical results for different values of 5^2 • The differences between 
the regulator-dependent values is, however, finite. Next, a straight-forward generalization of 
the discussion in Section fVI Dl tells us that in the absence of ^Vm^ing, the lattice regulator 
effectively fixes = c[ 2 /a s (here c[ 2 are the coefficients of the terms that multiply |tr P\ 2 
and |trP 2 | 2 in V(9) and that scale like l/a s at small a s ). Finally, in Section [Villi we show 
that for small a s , the Z^ symmetry is generically unbroken in our lattice calculation. This 
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means that c[ 2 > (otherwise there would be an instability). Therefore, if instead of letting 
MSDR fix bi t 2 = for us, we fix these coefficient such that the |trP| 2 and |trP 2 | 2 terms in 
V(6) have the same coefficients as they do on the lattice, we will find that the Zn symmetry 
is intact in MSDR as well. In particular, for small a s this means fixing b\ t 2 to have large 



and positive values. This of course is not surprising: Ref. H[ report a Zn — > Z2 symmetry 
breaking at m = and Nf = 1, but by increasing b\ and 62 to large positive values this 
surely will change and a Z N symmetric ground state will probably emerge. 

In the next section we depart from the EFT point of view, and simply study the lattice 
one loop potential as a function of its bare parameter space. Thus, we set the counter terms 
to zero and so fix 6 12 = 0. As a prelude to the full study of the phase diagram we first 
analyze the case of a s = a t in Section Wfll\ and indeed find that the Zn symmetry seems to 
be intact in the chiral limit. Next, in Section IVIII1 we fully explore the phase diagram of 
the potential Eq. (15.11) and find that the regime where the Zn symmetry is intact becomes 
extended when a s /a t is allowed to be different from one. Since in both cases, we find that 
the Zn is intact in the physically relevant regimes, we do not continue to ask what happens 
when we make SS c , t . nonzero. 



VII. THE PHASE STRUCTURE ALONG THE am AXIS: THE CASE OF SYM- 
METRIC LATTICE SPACINGS 

For symmetric lattice spacings we set X t = \ s and k s = K t = k into Eq. (15. ip . and turn to 
compare the values of the one loop potential for three vacua that realize the Zn symmetry 
differently. We do so for the physically interesting case of d = 4 (although below we present 
the analytic formulas for general d) and for different values of k (recall that at tree level 
massless fermions are obtained for k = l/2d - see Eq. (15.91) ). The vacua we considered are: 

• A vacuum denoted by that completely breaks the Z N symmetry. Here we set 

6 ab = 0, (7.1) 

and find 

ValN ^j (fY\J- ^ w \. (7.2) 
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A vacuum that preserves the Zjy symmetry. Here we set 



0" 



N 



(7.3) 



Substituting this into the one-loop potential and using 
and find 



(7.4) 



VzjN 2 



^] log 



S (d-2)/2 



-i 2V 



l-2Kd)+4K5j + 4/t 2 S 2 



> . 



(7.5) 



Note that here the integration is over a <i- dimensional Brillouin Zone and, correspond- 
ingly, S and S 2 are defined as sums over d terms: 

d d 



S = sin 2 k^/2, and S 2 = sin 2 k^. 



(7.6) 



Indeed, Eq. (17.41) is the way large-iV reduction embeds space-time into color space, 
and decompactifies the reduced direction. It is easy to check that Eq. (17.51) is exactly 
the one-loop potential one would obtain if one had an infinite lattice theory in all 
directions. 



We also studied a ground state with a Z 2 symmetry, i.e. that has 



e a 



ae[l,iV/2], 

Ti a E [N/2 + l,N\. 



(7.7) 



This means that out of the N 2 pairs of indices a and b there are iV 2 /2 that have a 
potential equal to Vq/N 2 , and the rest have an interaction given by 

r / j \ d-1 
I I fin \ 

AV Z2 



(S + l)(*-2)/2 



27r / I [((1 - 2nd) + Ak{S + l)) 2 + 4k 2 S 2 ] f 



(7.8) 



since for these 6 = it. Thus we see that this ground state has an energy of Vz 2 /N 
UV + AV Z2 )/N 2 . 
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An obvious uncertainty in our calculations is that we have only compared energies of the 
three vacua described above, and there may be other relevant vacua that we are ignoring. 

To obtain the phase diagram along the k axis, we scanned the values of Vz N ,Vz 2 , and 
Vq, in k G [0,2] and for Nf = 0.5, 1,2. 12 The integrations over k were done numerically 
with a trapezoid method whose grid had a resolution of 2ir/L in each direction with L = 
50, 100, 150. To obtain the L — > oo limit of these numerical integrations we performed linear 
extrapolations of Vz 2 , Vq and V Zn in (2%/L) p with p equal to the dimension of the Brillouin 
zone appearing in Eqs. (17. 2p . (17.81) . and (17. 5p (p = 3 in the former two and p = 4 in the 
latter). 

In Figs. (131), (HI), and (EJ) we present maps of the phase space along the k axis for the cases 
Nf = 0.5, 1, 2. For Nf = 1 we also zoom, in Fig. El on the regime of small k (that particular 
data set was generated for a single value of L = 80, but, in general, the variation with L 
was seen to be weak as long as we restrict to a s /a t = 1, as we do in this section). 



N f = 0.5 







■ None 

• Z N 
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FIG. 3: A map of the phase diagram symmetric lattice spacing and Nf = 1/2 (a single Majorana 
fermion in the continuum of the Zjq invariant phase) as a function of k (k = 1/8 is the chiral point 
and it has Zn symmetry intact). 

12 The case of Nf — 1/2 corresponds to a single Majorana fermion whose one-loop potential should vanish in 
the continuum limit of the chiral theory. The reason is simple: in that limit the theory is supersymmetric 
and the bosonic perturbative contribution must be canceled by the fermionic one. As is well known in 
the absence of the one-loop potential, non-perturbative instanton effects become important and make the 
ground state Z/v-symmetric. On the lattice, however, and away from the chiral limit, super-symmetry is 
broken, and we expect the one-loop potential to determine the ground state at sufficiently weak couplings. 
In that case the instanton effects should be exponentially small. 
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N f = 1 
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FIG. 4: As in Fig. but for N f = 1. 
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FIG. 5: As in Fig. U]but zooming on the regime of k G [0, 0.05]. 

As the figures show, close to k = one find Zn symmetry breakdown, which makes way 
to a Zn invariant state when we increase k. Surprisingly, this happens at quite small values 
of k: at k ~ 0.06 for N f = 0.5, at k ~ 0.04 for N f = 1, and at k ~ 0.01 - 0.04 at iV) = 2. 
In terms of the bare quark mass in lattice units these values correspond to am ~ 50 — 4. 
At even larger values of k, the Z^ symmetry breaks again. This, however, is less important 
since that regime corresponds to the so called 'super-critical' regime of Wilson fermions 
which one needs to avoid in lattice simulations (it is not in the same universality class of 



QCD [2J). 



Importantly, the chiral point (and a generous vicinity thereof) is at k — 1/8 and this is 
within the Z^ symmetric phase for all choices of Nf. Also, as anticipated, the Z^ symmetric 
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N f = 2 
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FIG. 6: As in Fig. but for N f = 2. 
phase becomes more extended with increasing number of flavors. 

VIII. THE PHASE DIAGRAM IN THE PLANE OF a t m AND a s /a t . 

We now turn to map the phase diagram in the plane of a t m and a s /s t . As in Section IVHl 
we performed the numerical integrations with a fixed grid in momentum space (which in 
this section was set to be 2n/L with L = 60, 90) and extrapolated to zero grid spacing (see 
discussion in previous section). In contrast to the case of a s /a t = 1 discussed above, when 
a s /at is small this extrapolation is important to perform and sticking to a fixed value of L can 
result in an erroneous phase diagram. For example, in Figs. (EJ) — (JHJ) we plot the difference 
in the potentials (Vz N — Vz 2 )/N 2 for a t m = 0, a t /a s = 0.01 and Nf = 0.5,1,2, versus 
(27r/L) 3 . While Vz N and Vz 2 are expected to depend linearly on {2n/L) A and (27r/L) 3 , 
respectively, then the fact that the plots are linear means that most of the variation in 
the difference Vz N — Vz 2 reflects the variation of Vz 2 . Our data was generated for L = 
70, 90, 100, 140, 160, and as is clear from the figures, at L — > oo we see that Vz N < Vz 2 and so 
that the Zn symmetry is intact. Fortunately, we see that the linear behavior sets in already 
at the small values of L where we can perform the numerical integration at a reasonable 
computational cost. Therefore, from here on we shall restrict ourself to performing the 
numerical integrations with L = 60 and 90, and map the phase diagram according to the 
L — > oo linear extrapolations of these potentials. Let us emphasize, however, that while it is 
important to perform the large-L extrapolations at small values of at/a s , then at moderate 
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FIG. 7: The difference in energy between the Zj^ invariant state and the Zi invariant state for 
Nj = 1/2, atm = and a s /a>t = 0.01 versus the momentum space resolution used to perform 
the numerical integrations over the Brillouin zone. At L — > oo the difference is negative and the 
ivTv- invariant vacuum is energetically preferable. 

values of this parameter the results we obtain prior to the extrapolations are quite close to 
their large- L limit. This is expected since at very small values of a t /a s , only the vicinity 
of the Brillouin Zone origin is important, and a finer grid is necessary. In practice, we find 
that it is only at a s /a t < 0.2 that the linear extrapolations are important, while for larger 
a s /a t , a numerical integration with L = 90 is already reflecting the situation at L = oo. 

The results we find are quite interesting: we see that introducing an anisotropy makes 
the range in which the Zn symmetry is intact more extended. For physically relevant values 
of the quark mass, however, nothing dramatic happens and the ground state is still Zn 
invariant. We present the map of the phase diagrams in Figs. (fT0l) - (jT2l) . 

IX. CONCLUSIONS 



In this paper we analyzed the large- N volume independence of four dimensional QCD 
with adjoint fermions, and we find that it works in weak coupling if one regularizes the 
theory on the lattice with Wilson fermions. Specifically, we studied this regularization for 
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FIG. 8: The same as in Fig. but for N f = 1. 
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FIG. 9: The same as in Fig. El but for N f = 2. 



different number of flavors, different quark masses, and a varying anisotropy between the 
spatial and temporal lattice spacings. Our calculation is performed at one-loop and we 
calculate the corresponding effective potential V as a function of the eigenvalues {e t6a } of 
the holonomy in the reduced direction (the Polyakov loop P). We find that V{9) prefers 
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FIG. 10: A map of the phase space according to the one loop potential for anisotropic lattice 
spacings and Nf = 1/2 (a single Majorana fermion in the continuum of the Zpj invariant case) 
as a function of a^m and a s /at- Bursts (red) denote the points in the parameter space where the 
Ztv symmetry was seen to be completely broken, circles (green) denote the regime where the Zn 
symmetry is broken down to Z2, and dots (blue) the regime where the Zn symmetry is intact. 
Note that the scan in the proximity of atva = was done with a finer resolution in a^m. 



a Zn invariant ground state for moderately light (and even quite heavy) fermions, and for 
most values of the anisotropy. Our results hold for 1/2, 1, 2 Dirac flavors (one, two and four 
Majorana fermions in the continuum). If, however, the fermions are extremely heavy, then 
the reduced theory spontaneously breaks its Zn symmetry and large- TV volume independence 
breaks down. We also see that the Z^ symmetry is broken in the (physically uninteresting) 
super-critical regime of the lattice Wilson theory. 

Another goal of this paper was to understand whether there is any tension between the 
results of Ref. which treated the volume-reduced model as a 3D continuum effective 
field theory (EFT), and those of Ref. 2j (that worked directly with the AD gauge theory 
in the continuum), and whether the results of the former signal any problems with large- N 
reduction. In particular, for massless quarks Ref. Q finds that V(9) has a Zn invariant 
vacuum, while Ref. |ll] sees a vacuum that breaks the Z^ symmetry down to Z 2 - In fact, 
from the previous paragraph it seems that the results in Ref. llj also contradict the results 
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FIG. 11: As in Fig. [TO] but for N f = 1 
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FIG. 12: As in Fig. [TU] but for Nf = 2. 



of our lattice calculation, in which the 3D spatial dimensions of the reduced model are 
discretized to have a finite lattice spacing a s . 

Before we describe how we resolved these 'contradictions', we wish to emphasize that 
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treating the reduced model as a three-dimensional EFT may be useful but is not necessary. 13 
In particular, in the lattice calculation performed in this paper we treat the L t — 1 reduced 
model as a theory that is defined with a fixed cutoff, or equivalently that has finite bare 
lattice parameters. For values of lattice parameters where the center symmetry is intact, the 
theory is large- N equivalent to a corresponding theory with the same bare lattice parameters, 
the same field content, and the same cutoff, but with L t = oo. This equivalence is true not 
only at low energies, but all the way up to the cutoff scale. Removing the cutoff a s from 
this construction is finally done after taking the large-N limit, and according to the RG 
equations of the 4D large- N theory. 

While our approach is the standard way one defines large- N reduction nonperturbatively 
[26], it does not mean that the approach of Ref. jlj| is not useful, and we still need to 
understand how to resolve the apparent contradiction between the results of this approach 
and what we and Ref. [2j find for the symmetry of the ground state. For that purpose we 
first showed in Section IVII that the three-dimensional EFT defined by the reduced model is 
non-renormalizable. One of the consequences of this is a set of linear divergences in V(9) 
that depend on 9. These are radiatively generated as mass terms for the Polyakov loops 
and from the point of view of EFT we need to cancel them. Thus, one needs to add to 
the action of the EFT certain counter terms. These turn out to be relevant operators and 
after canceling the UV divergences, they leave us with extra finite additions to the one loop 
potential that depend on new low energy constants (LEC). Specifically, we showed that the 
finite contribution of the counter terms is of the form 

^finite- J a*x{b x \X,vP{x)\ 2 + b 2 \X,vP 2 (x)\ 2 }. (9.1) 

where the LEC b\ and 62 have mass dimension three. We also identified the counter terms 
for general values of L t and of the spatial dimension D (see Section IVTl) . Interestingly, the 
counter terms constitute a subset of the operators that Ref. [8| suggested to add to the 
Eguchi-Kawai model in order to make its ground state Zjv-invariant. In this paper we see 
that in the EFT approach, their presence is dictated by the regularization process of the 
reduced model. (Note that if the theory has a ground state with a Zjsi symmetry, then these 



13 



The reason why Ref. ll| choose to do so is because it may open a window for using three dimensional 



analytic techniques to study the four dimensional theory. 
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operators, when projected back to AD, do not affect the leading large- N dynamics - again 
see Ref. [3]). 

Therefore, the reduced model, seen as a three-dimensional EFT, is defined not only by 
parameters like the quark mass m and the number of flavors Nf, but also by 61 and 62- This 
tells us what is the cause for the different results one obtains in different regularizations: 
they reflect different choices for b\ and b 2 . To make this clear let us detail the choices made 
for these LEC by the different calculations. 



I. The three dimensional EFT continuum calculation of Ref. H[ : 



Here the subtraction of the infinities was done with minimally subtracted dimensional 
regularization (MSDR). This choice effectively sets power law divergences to zero and 



replaces them by finite 0-dependent functions, 
was not added to the action studied in Ref. 



he fact that 5S^ n i te from Eq. 09. ip 
111 ] means that MSDR actually fixed 



bi = &2 = 0. This choice is made implicitly by the regulator and is why MSDR is 



considered 'dangerous' when dealing with power-law divergences 23l |. 

II. The lattice calculation presented in our work (Section IV]) : 

Above we emphasized we do not treat the reduced model as an EFT. Despite this, and 
to understand the difference between our results and those of Ref. Ill ], we momentarily 
choose to depart this approach, and think about our lattice as a regulator for a 3D 
EFT which is alternative to MSDR. It is easy to see that, from this point of view, the 
lattice fixed b\ and 62 to specific values that scale like l/a s for small a s (because we 
did not add Eq. (19. II) to our lattice action - for details see Section IVIDI) . 

There is absolutely no reason to expect that the implicit choices made for b\^ i n MSDR 
and on the lattice should lead to the same physical results. In particular, it is quite possible 
that these choices lead to different realizations of the symmetry in the ground state. 
Indeed this is what happens: the choice made by MSDR tends to break the symmetry, 
while the one of the lattice (with Wilson fermions) tends to preserve it. 

To obtain the same physical results in different regularizations, one would need to fix 
the physical parameters (m, Nf, etc.), explicitly introduce the terms in Eq. (19. ip into V(8), 
and tune the values of bi and 62 i n a regularization dependent way (instead of letting the 
regulator fix them implicitly). Then, there will be some choice within MSDR, 6 12 = 6^| DR , 
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that will give the same physical results that a different choice, b\$, = fr 1 ^ 106 , gives on the 
lattice. 

Moreover, to make a regularization-independent statement on the absence of Zn symme- 
try in the ground state of the theory, one needs to show that, within a certain regularization, 
the center symmetry breaks for all values of b\ and hi. As we say above, this was not in- 



cluded in the analysis of Ref. 11( and so the result of that work does not mean that large- N 



reduction of QCD with light adjoint fermions is invalid. Put differently, Ref. [ll| effectively 
studied V{6) in a subspace of the full parameter space of the EFT; a subspace that turns 
out to have a broken Zn symmetry for physically relevant values of the quark mass. If we 
view our lattice calculation as an EFT, then the choices made for 6^2 by the lattice can also 
be seen as a restriction to a subspace in parameter space. But in contrast to what happens 
with MSDR, in our case the subspace defined by 5\ a 2 tlce turned out to generically have a 
^Tv-invariant ground state for V{6). 



Is it possible to change the results of Ref. rjj by making its choice of 61,2 nonzero? The 



answer seems to be yes. For example, for m = and Nf = 1, Ref. H] reported a breakdown 
of Zjy — > Z2, and it is fairly clear that by increasing both b\ and 62 to sufficiently positive 
values, one can get rid of this symmetry breakdown. Indeed, this fact is what makes the 
results of Ref. consistent with the ones we present in the current paper, namely that 
large-iV reduction works for a physically relevant range of the parameters in the gauge 
theory. This also means that there is no tension between the results of Ref. [ll] and those 
of Ref. 0. 

The fact that we can view the lattice definition of the L t — 1 reduced model as a regular- 
ization of a three-dimensional nonrenormalizble EFT with fixed b\p teaches us the following 
important lesson. Other lattice constructions of the L t = 1 model (for example ones similar 



to [25] which use staggered fermions, or any other type of fermions) can be seen as alter- 
native regulators for the same EFT. Then, if we do not add the terms in Eq. (19.1 p to their 
action, they will implicitly choose their own values for &i,2- Importantly, it is not guaranteed 
that these choices will generically lead to a Z^ symmetric ground state in the physically 
interesting regime. In that sense, the result we present in this paper for Wilson fermions 
cannot be anticipated in advance. Clearly, this means that it will be wise to perform one- 
loop calculations of the form we did in this paper for each regularization of a single-site 
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model, prior to its (computationally costly) numerical Monte-Carlo study. 



The question of whether the symmetry is intact in our regularization also at mod- 
erate couplings, where the one-loop calculation is unreliable, can be answered only via 
non-perturbative Monte-Carlo simulations. For example, in the Nf = 2 case, the study in 
Ref. [25J suggests that the answer is rather complicated and may be sensitive to the bare 
lattice parameters. The results we present in Section IVIIII on the phase structure in the 
a s /at — dtm plane can be viewed as another example of a moderate case of this sensitivity, 
which is harmless for the large- N reduction program. 



Indeed, in the companion publication [l2j non-perturbative Monte-Carlo simulations were 
used to explore large- N reduction and find evidence that for some values of the lattice 
coupling and quark masses, the theory with symmetric lattice spacings can be successfully 
reduced to a single site in all the euclidean directions. Anticipating this result using a 
one loop calculation is, unfortunately, not straightforward, and compared to the one-loop 
calculation presented in this paper, is complicated by IR divergences. These need to be taken 
into account in a similar fashion to the way Ref. [23] estimated the free energy of singular 



tolerons. Nonetheless, the values of the quark mass at which we see a transition from a Z 



N 



broken phase into a Zn invariant phase in the Monte-Carlo simulations of Ref. 12j are in 
qualitative agreement with the values we find analytically in the current paper. 
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